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We present a both ultraviolet and infrared regularization independent analysis in a symmetry pre- 
serving framework for the N = 1 Super Yang-Mills beta function to two loop order. We show 
explicitly that off-shell infrared divergences as well as the overall two loop ultraviolet divergence 
cancel out whilst the beta function receives contributions of infrared modes. 

PACS numbers: ll.10.Gh, 11.15.Bt, 11.30.Pb 



CD 



> 

On 
(N 

G\ 
O 

o 



X 



I. INTRODUCTION 



In this contribution we aim at shedding some 
light, from the perturbative standpoint, on the ori- 
gin of loop contributions to the beta function of 
N = 1 Supersymmetric Yang-Mills (SYM) theo- 
ries. We employ an automatically invariant, four 
dimensional (minimal) framework in which neither 
Ward identities need to be imposed as constraint 
equations nor modifications in the original physi- 
cal Lagrangian are performed. Furthermore both 
finite and divergent parts of the Feynman ampli- 
tudes are displayed which is particularly useful to 
discuss the subtleties involved in the corrections 
of the beta function beyond one loop order in this 
model as we shall discuss later on. 

The beta function of SYM theory has been a 
subject of controversy in both non-perturbative 
and perturbative calculations. One of the debates 
has been named "anomaly puzzle". Because the 
U(l)n current is in the same multiplet as the trace 
of the energy momentum tensor, chiral anomaly 
and trace anomaly are equally subjected to the 
Adler-Bardeen theorem [l| . As the beta function is 
proportional to the trace of the energy momentum 
tensor, one could assert on general grounds that 
contributions beyond one loop order were zero. On 
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the other hand Novikov, Shifman, Vainshtein and 
Zakharov (NSVZ) were the first to obtain an ex- 
act expression for the beta function Q. In addi- 
tion, perturbative calculations using dimensional 
reduction |3| revealed higher loop contributions to 
the beta function. Other regularization methods 
have obtained corrections beyond one loop as well 
01 , @ i ■ Arkani-Hamed and Murayama Q have 
shed some light on the anomaly puzzle providing a 
simple and elegant explanation. They have shown 
that the vector multiplet does not have canonical 
kinetic terms after dilations (trace anomaly) which 
can only appear after an additional rescaling is per- 
formed. Hence the anomaly from this modified di- 
lation will no longer belong to the same multiplet 
as the U(1)r anomaly. As a result, the beta func- 
tion corrections are not constrained by the Adler 
Bardeen theorem. However, some controversies re- 
main Q. 

The other controversy on this matter con- 
cerns the possible infrared mode contributions to 
the beta function which again affects both non- 
perturbative and perturbative calculations. NSVZ 
[2| derived the exact beta function in a frame- 
work based on instanton analysis which borne out 
its infrared origin. Contrariwise, Arkani-Hamed 
and Murayama Q claim, within a purely Wilso- 
nian framework freed from infrared subtleties, that 
the exact beta function depends only on ultravio- 
let properties of the theory. This discussion also 
appears in a perturbative analysis. Within the 
supergraph approach to supersymmetric models, 
along with on-shell infrared divergences of Yang- 
Mills theory, additional off-shell infrared diver- 
gences appear Q which must be consistently sepa- 
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rated from ultraviolet ones before renormalization 
is effected. The unclear distinction between these 
two types of infinities is at the heart of this de- 
bate. The answer to this matter passes through 
an unambiguous distinction between the infinities 
involved and the arbitrary scales which are byprod- 
ucts of the subtractions. In dimensional reduction 
9], [10] the two loop correction to the beta func- 
tion appears from a local evanescent operator typ- 
ical of the method. Such operator is absent should 
the calculation stay in the physical space-time di- 
mension. Therefore Grisaru, Milewski and Zanon 
themselves conjectured that no divergence should 
occur beyond one loop. This is correct as we ver- 
ify later in this contribution. However it does not 
mean that the two loop beta function vanishes but 
instead that the naive perturbative derivation of 
the beta function based on renormalization con- 
stants needs some reinterpretation and modifica- 
tion. Such modification appears to be related to 
scaling anomalies [ll| . It is important to observe 
that in dimensional regularizations the treatment 
of amplitudes which are simultaneously infrared 
and ultraviolet divergent is involved because they 
can be mixed. Mas, Perez- Victoria and Seijas have 
evaluated the two loop contribution within Differ- 
ential Renormalization [4] which does not recourse 
to a dimensional extension on the space-time di- 
mension. In this framework no infinities appear 
by construction. They claim that the beta func- 
tion two loop correction depends on infrared di- 
vergences yet the latter is seen as playing a pas- 
sive role. We believe that a calculation in which 1) 
no dimensional continuation on space-time is per- 
formed; 2) infrared and ultraviolet divergences de- 
couple clearly; 3) one distinguishes the correspond- 
ing scales without the need of introducing extra 
parameters to ensure their independence; 4) gauge 
and supersymmetry invariance are maintained at 
every loop order; 5) only one renormalization scale 
appears as divergences are subtracted and 6) one 
can keep track of both finite and infinite parts of 
the amplitude, can shed more light on the pertur- 
bative analysis of this problem. 

Implicit Regularization is a good candidate for 
such a task. It has been proved to be invari- 
ant to all orders if well known surface terms are 
set to zero which is closely related to momen- 
tum rou ting invariance in the Feynman diagrams 
E3.E3.Il3. El- Ultraviolet infinities need not be 
evaluated and remain as a systematic set of basic 
divergent integrals which depend upon a unique 
renormalization scale A. Their derivatives with re- 
spect to that scale are also basic divergent inte- 



grals. Off-shell infrared infinities can be likewise 
defined save that the basic divergent integrals de- 
pend on an infrared arbitrary scale which is inde- 
pendent of A. Moreover, it was shown that Im- 
plicit Regularization complies with locality, uni- 
tarity and Lorentz invariance because the Bogoli- 
ubov's recursion relations which locally executes 
the forest formula to subtract general subdiver- 
gences can be implemented in such a scheme. 

In the present contribution we calculate the two 
loop beta function of N — 1 Super Yang-Mills the- 
ory in an ultraviolet and off-shell infrared regular- 
ization independent way in the framework of Im- 
plicit Regularization. We verified in our calcula- 
tion at two loop order the conjecture by Grisaru, 
Milewski and Zanon E3 that after subtraction of 
subdivergences no ultraviolet infinities occur be- 
yond one loop in N=l Super Yang-Mills when 
space-time dimension is not modified. More im- 
portantly we conclude that the two first terms in 
the loop expansion of the beta function receive con- 
tributions related to on-shell infrared divergences. 



II. PURE N=l SUPERSYMMETRIC 
YANG-MILLS 

Consider the standard pure N = 1 supersym- 
metric Yang-Mills classical action E3 

Ssym = J d A xd 2 e TrIF 2 . 

We work with the supersymmetric background 
field method where gauge invariance is kept man- 
ifest. The one loop contributions depicted by the 
diagrams of figure 1 read, after performing super- 
symmetric algebra, 

d a = -\c A J d i pd i ew a ( P )r a (-p) 

X J (2tt) 4 PCp + fc) 2 

D B = i\c A |d 4 pd 2 (?W (p)W a (-p)p 2 

/■ A - d 4 k 1 
X J (2^)4fc 4 (fc-p) 2 

where £ is a gauge parameter and the tilde in- 
dicates the Fourier transform. The superscripts 
A U v , Ai r indicate that the integrals are implicitly 
regularized in the ultraviolet and infrared limits 
so to allow for algebraic manipulations with the 
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FIG. 1: One loop diagrams of the two-point function. 



integrand. Hereafter, for brevity, we adopt the no- 
tation J fc = / A "" ,A4r ^. 

By simple power counting we conclude that the 
first integral in k is only ultraviolet divergent while 
the second one is only infrared divergent. In order 
to avoid as long as possible contact with infrared 
regularization schemes we leave the infrared inte- 
grals in their original form assuming an adequate 
regulator that allows for algebraic manipulations. 
In the appendix we specify two concrete infrared 
regularizations. We can write Da in terms of the 
classical action after integrating over half of the 
supercoordinates, to obtain 



Da = -^C A i J d 4 p d 2 9 W a (p)W a (-p)I(p) (1) 
D B = i\c A ! d A p d 2 9 W a ( P )W a (-p)p 2 U(p) 



(2) 



where we have defined for later use 



kf 



The ultraviolet divergent integral I(p) can be 
written in the convenient form using Implicit Reg- 
ularization (see for instance ITU S 



J(p) - I log (y)-b\n[-^ 



2b 



Ii og (X 2 ) + lf in {p 2 ,X 2 ) 



(3) 
(4) 



where b 



(Itt) 



A is a nonvanishing arbitrary pa- 



rameter which plays the role of renormalization 
group scale and 



Ilog{>?) = 



1 



(k 2 - A 2 ) 2 



(5) 



is a one loop basic divergent integral [20j. 

The two loop contributions depicted in figure 2 
can be written, omitting a common factor — ^ — -, 




FIG. 2: Two loop diagrams of the two-point function. 



as 



D 1 = / d 4 pd 4 9W (p,9)W a (-p,9)a» 



Do 



k,q k 4 (p~k) 2 q 2 (k-q) 2 
d 4 9 d 4 p r"{-p, 0)f"(p, 9) 



Pp,Pv 



k,a q 4 {p + q) 2 k 2 {k + qf 



+4 



(p + q)v<in 



k,q q 4 (p + q) 2 k 2 ( g + k) 2 



Da 



- J d 4 9d 4 p T»{~p,9)T»{p,9) 

(p + q)n(p + k) v 



k q q 2 k 2 (p + q) 2 (p + k) 2 (k — q) 2 



(p + q) tJ .(p + k)u 



kq {P + q + k) 2 (p + k) 2 q 2 k 2 (p + q) 2 



Da 



d 4 9 d 4 p r"(-p, 0)T v (p, 



k,q (q + kYq 2 k 2 - 

Di, D% and D4 are both ultraviolet and infrared 
power counting divergent while D3 is only ultra- 
violet divergent. In addition D±, D2 and -D4 con- 
tain ultraviolet subdivergences which must be sub- 
tracted. 

We start by evaluating Di 



Di 



d 4 p d 4 9 W {p, 9)W a {-p, 9)a 



Oft 



(P - k)» 
k k 4 {p-kf 



I(k). 



The ultraviolet subdivergence is removed mini- 
mally by subtracting Ii og {\ 2 ) from I(k). We de- 
fine D[ by substituting I(k) -> (p 2 , A 2 ) in the 
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expression above. Then, after using the identity 
V Q W Q = V a W a , D[ becomes 



D' 



d 4 P d 2 e w a ( P ,8)w a (- P , e)p 2 



i t k'(p-kr- 

It is straightforward to see that D 2 can be cast 



as 



ih = i d 4 e d* P r»(- P ,d)r»(p,e) 



P^Pv 



+4 



q q 4 {p + q) 2 
q v q^ 



±Pv 



q^(q) 

\ q 4 (p + qf 



q 4 (p + q) : 



ri(q) 



and again D' 2 is obtained from D 2 substituting I(q) 
by P in (g 2 ,A 2 ). 

Consider now the amplitude D3 which contain 
no subdivergences. After some simple algebra we 
find 

d' 3 = d 3 = IJ d 4 e d 4 P f»{- P ,6)f»(jp,e) 

V -^J{p) + 2I, v 



where 



kuQv 



k J q k 2 q 2 {k-pf (k-qf {q - pf 
and the finite integral [l 

J(P) 



k J q k 2 q 2 (k - pf (k - qf (q - p) 2 
= 6C(3) 
(4tt) V ' 

As for D4 after changing variables q — » q 1 = k+q 
we separate it in ultraviolet and infrared divergent 
pieces by multiplying and dividing by (p + q) 2 . It 
results 

At = - 9 -f [d 4 6d 4 pr»(-p,9)T»(p,d) 



p 

2p a 



Hq) 



q 4 (p + q) 2 J q q 2 {p + q) 2 

q a I(q) 



q 4 {p + qf 



D' 4 is obtained as before replacing I(q) by 
lf m (q 2 ,X 2 ). 



The following integrals are contained in the ex- 
pressions for -D-'s. We present the results [14{ up to 
surface terms which we systematically set to zero 
to implement gauge (and momentum route) invari- 
ance as discussed in the introduction fl4| : 
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q 4 (q + 
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q 4 (q + 


P) 2 


q^q 




q 4 (q + 


pf 


b p^lPv 
2 p 2 ' 




Qu.Q 


J 



In - 



A 2 



/S(A 2 ) + 61n(- 



A 2 



-b P -\n(^ 



A 2 



9^ 
4 



/ ; „ 5 (A*)-61n(-^ 



2b 



• 9 



2 
P 



4 

2\ b 



4 2 i(A 2 ) 



q 4 (q + pf 



6 6 / p 2 
4 + 2 ln 



V" - 4 



^b 2 -b 2 ^ 
3 9 



bhog{X 
2 



b 2 In 



El 

A 2 



-J(P) 



r p2 



where / ; (2 j(A 2 ) = J feW ^l n (-%^] and 
Iiogi^ 2 ) is defined in JS]). 
The partial sum 



D' 2 + D 3 + D' 4 = ± I d l Hd l i,f"i-,,.0)f"{i>. 
x(^-pV)| ; , y)J 



^ m (<Z 2 ,A 2 ) 6 2 C(3) 



26 / , . / p 2 \ 8 . 6tt 
-61n( I + -6- 



3G 



(0) 



is transverse as it should because of gauge invari- 
ance which is manifest in the background field 
method, showing that our calculation is consistent 
(invariant). Therefore it can be written in terms 
of the classical action by integrating over half of 
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the supercoordinates to yield 



D' 2 + D' 3 + D' 4 = -j d 2 9 d 4 p W a (-p, 9 
P' m (g\\ 2 ) b 2 ((3) 

q q 4 (p + q) 2 p 2 
2 ^ 8. bn 2 ' 




x 2 



36 



(7) 



Now we can add all the contributions to the two 
loop two-point function as 

4 r ~ ~ 

J2 D > = J d2 ° d4 p w a (-p,o)w a (p,e) 



3g 2 C 



A 



P in (fc 2 ,A 2 ) 



In 



V J k k 4 (p-k) 2 
A 2 



+ i b C(3) - 2b 2 + 



24 



.(8) 



Notice that the ultraviolet divergences at two 
loop order cancel out magically in the sum D[ + 
D' 2 + D' 3 + IX as predicted by Grisaru, Milewski 
and Zanon [1QJ. As we have said in the introduc- 
tion, this does not mean that the two loop beta 
function coefficient vanishes. In order to evaluate 
the beta function one can use the renormalization 
group equation. Adding the tree, one loop and two 
loop contributions given by ([1]) , ^ and © allows 
us to write the (ultraviolet) renormalized two-point 
Green function as 



G (2) („2N = J_ 



3g 2 C 2 



36 In (-^J + 6b + & 2 U( P ) 
2bp 2 U(p) + bp 2 U (2) (p) 



36 2 C(3) 



6 2 ^ 2 



G 



(9) 



where the first term is the tree level contribution, 
U(p) was defined in the one loop calculation and 
U^(jp) stands for 



U (2 \ P ) 



k k*(k- P y n v a 2 



(10) 



Notice that U^(p) is off-shell infrared divergent 
and it contains both an off-shell infrared (A) and 
ultraviolet (A) scales (see the appendix). They are 



obviously independent by construction. Moreover 
one can easily check that 



^(Gi 2 J n (P 2 



which guarantees the independence of the scales 
A and A. Here R is the operation that subtracts 
off-shell infrared infinities. We leave the off-shell 
infrared divergences (implicitly regularized) repre- 
sented by U(p) and U^(p) in an integral form to 
show explicitly that they play no role in the eval- 
uation of the beta function. 

G^renip 2 ) obeys the renormalization group equa- 
tion [21 : 



= 0. 



5=0 



(11) 

j^g 2 + 0(g 4 ) [J and © we can 



Using 7 4 = 

solve the above equation for /3(g). Writing (3(g) 
hg 3 + b 2 g 5 + 0(g 7 ) we find 



bi 



3 C A 



4 (4tt) 2 



(12) 



b 2 = tA 



bp 2 X 2 ^U^(p)-2b 2 



bp 2 U(p) 



(13) 



Note that the second term on the right hand 
side of (|13[) comes from the term proportional to 
ln(-p 2 /A 2 ) of the last line of ©. The terms 
U^ 2 \p) and U(p) are separately off-shell infrared 
divergent. However \ 2 -^U^ 2 \p) = —U(p) inde- 
pendently of the infrared regulator and thus 



b 2 = — 



r 2 



4 (4tt) 



(14) 



It is important to note that such cancelation does 
not depend on the infrared regularization and since 
just U^ 2 \p) and U(p) are (off-shell) infrared diver- 
gent the two first universal coefficients of the beta 
function are independent of off-shell infrared in- 
finities. Indeed in the appendix we explicitly eval- 
uate the integrals U(p) and U^ 2 \p) in two different 
infrared regularizations showing such cancelation. 
The above beta function coefficients agree with the 
literature (see for instance [H). 



III. CONCLUDING REMARKS 

We used, we believe for the first time, a four 
dimensional invariant framework in which gauge 



G 



and supersymmetry invariance are automatically 
maintained, and both finite and infinite parts of 
the amplitudes are made explicit throughout the 
calculation, without any modification to the origi- 
nal effective action (i.e. with the simplest Feynman 
rules), to study the loop contributions to the beta 
function of N = 1 SYM theory. In the Feynman 
gauge, along with on-shell infrared divergences, ad- 
ditional off-shell divergences appear as expected 
Q; they must be carefully dealt with, in a way 
that they do not get mixed with ultraviolet ones. 
This is achieved through the identification of dis- 
tinct scales in connection with the R and/or R 
operations. 

As conjectured in fiol ] for four dimensional in- 
variant methods, we verified within our framework 
that after subdivergences are subtracted, no overall 
ultraviolet divergences survive at two loop order. 
However the beta function does not vanish as cal- 
culated with the finite part of the amplitude which 
obeys a regular renormalization group equation. 

As it becomes clear from our equation (|13l) . off- 
shell infrared divergences cancel out automatically 
in a (off-shell infrared) regularization independent 
fashion (that is to say, no R operation is neces- 
sary) and hence there is no contribution what- 
soever stemming from them. We have explicitly 
verified within two off-shell infrared regularization 
frameworks the cancelation of the first and third 
term in (|13l) using in one case a dual version of Im- 
plicit Regularization in coordinate space and in the 
other case a simple fictitious mass in the propaga- 
tors (see appendix) . However the term in the mid- 
dle, proportional to b 2 , stems from the derivative 
of a non local on-shell infrared divergence propor- 
tional to ln(— p 2 /X 2 ) which multiplies the W(p, 9) 
fields and is integrated over the whole p four di- 
mensional space. Therefore it becomes apparent 
that (on-shell) infrared modes contribute to the 
two loop beta function. This is somewhat not sur- 
prising as it may happen in other supersymmetric 
theories [l9j |. 

It is clear that a naive perturbative derivation of 
the beta function from renormalization constants 
must be modified taking into account thepresence 
of the rescaling anomalies discussed in [Tl|. Our 
approach reveals itself as very adequate for such 
analysis. 
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APPENDIX 

In [rj| the integrals U (p) and (p) were evalu- 
ated using a dual version of Implicit Regularization 
that deals with infrared divergences. Using their 
results: 



U = 



i ^ 09 (A- 2 ) + Mn^-g | + 2b 



c/( 2 




1 T {2) (\- 2 \ 



where A is an infrared scale independent of the 
ultraviolet one, A = (4/e 27 ) A with 7 being the 
Euler-Mascheroni constant. Here Ii og (\~ 2 ) and 
/^(A -2 ) are infrared basic divergent integrals 
(analogous to the ultraviolet ones). Substituting 
these results in ([T^]) we see the cancelation of the 
infrared divergent pieces and the nonlocal terms 
and we get the correct value of 62 (HU- 

Alternatively one can introduce an infrared mass 
regulator /1 in the propagators finding: 



U = 



U 



(2) 




Again substituting these results in (IT31 we get the 
correct result of 62- 
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